DARK MATTER PARTICLES 



V.Berezinsky 

INFN, Laboratori Nazionali del Gran Sasso, 67010 Assergi (AQ), Italy 
and Institute for Nuclear Research, Moscow 

The baryonic and cold dark matter are reviewed in the context of cosmological 
models. The theoretical search for the particle candidates is limited by supersym- 
metric extension of the Standard Model. Generically in such models there are just 
two candidates associated with each other: generalized neutralino, which compo- 
nents are usual neutralino and axino, and axion which is a partner of axino in 
supermultiplet. The status of these particles as DM candidates is described. 



1 Introduction 

Presence of dark matter (DM) in the Universe is reliably established. DM in 
the form of compact microlensing objects (machos) are directly observed in the 
halo of our Galaxy by MACHO, EROS and OGLE collaborations. Rotation 
curves in our Galaxy and in many other galaxies provide evidence for large 
halos filled by nonluminous matter. The virial (gravitational) mass of clusters 
of galaxies is about ten times larger than their luminous masses. IRAS and 
POTENT demonstrate the presence of DM on the largest scale in the Universe. 

The matter density in the Universe p is usually parametrized in terms of 
fl = p/pc, where p c « 1.88 • 10~ 29 h 2 g/cm 3 is the critical density and h is the 
dimensionless Hubble constant defined as h = H / (100km. s^ 1 .Mpc^ 1 ). Dif- 
ferent measurements suggest generally 0.4 < h < 1. The recent measurements 
of extragalactic Cepheids in Virgo and Coma clusters narrowed this interval 
to 0.6 < h < 0.9. However, one should be cautious about the accuracy of this 
interval due to uncertainties involved in these difficult measuremets. 

Dark Matter can be subdivided in baryonic DM, hot DM (HDM) and cold 
DM (CDM). 

The density of baryonic matter found from nucleosynthesis is given El as 
0.009 < ilbh 2 < 0.02. The baryonic cosmological density provided by the mass 
of intracluster gas is very close to this value, fl^ Lust h 3 ^ 2 « 0.050 and references 
therein. 

The structure formation in Universe put strong restrictions to the proper- 
ties of DM in Universe. Universe with HDM plus baryonic DM has a wrong 
prediction for the spectrum of fluctuations as compared with measurements of 
COBE, IRAS and CfA. CDM plus baryonic matter can explain the spectrum 
of fluctuations if total density fio ~ 0.3. 

There is one more form of energy density in the Universe, namely the 



1 



vacuum energy described by the cosmological constant A. The corresponding 
energy density is given by £Ia — A/(3Hq). QuasarJensing restricts the vacuum 
energy density: in terms of 17a it is less than 0.7u . 

Contribution of galactic halos to the total density is estimated as 51 ~ 
0.03 — 0.1 and clusters give f2 ps 0.3. Inspired mostly by theoretical motivation 
(horizon problem, flatness problem and the beauty of the inflationary scenar- 
ios) O = 1 is usually assumed. This value is supported by IRAS data and 
POTENT analysis. No observational data significantly contradict this value. 

There are several cosmological models based on the four types of DM de- 
scribed above (baryonic DM, HDM, CDM and vacuum energy). These models 
predict different spectra of fluctuations to be compared with data of COBE, 
IRAS, CfA etc. They also produce different effects for cluster-cluster corre- 
lations, velocity dispersion etc. The simplest and most attractive model for 
a correct description of all these phenomena is the so-called mixed model or 
cold- hot dark matter model (CHDM). This model is characterized by following 
parameters: 

ClA = 0, Oo = fi& + QcDM + &HDM = 1, 

H PS 50 kms^Mpc^ih ps 0.5), 
SIcdm ■ Vhdm ■ n b « 0.75 : 0.20 : 0.05, (1) 

Thus in the CHDM model the central value for the CDM density is given by 
QcDNih 2 = 0.19, with uncertaities within 0.1. 

The best candidate for the HDM particle is r-neutrino. In the CHDM 
model with n v = 0.2 mass of r neutrino is m„ T ~ 4.7 eV. This component 
will not be discussed further. 

The most plausible candidate for the CDM particle is probably the neu- 
tralino (x): it 1S massive, stable (when the neutralino is the lightest super- 
symmeric particle and if R-parity is conserved) and the xx-armihilation cross- 
section results in Q x h 2 ~ 0.2 in large areas of the neutralino parameter space. 

In the light of recent measurements of the Hubble constant the CHDM 
model faces the age problem. The lower limit on the age of Universe to > 13 Gyr 
(age of globular clusters) imposes the upper limit on the Hubble constant in 
the CHDM model H < 50 kms^Mpc' 1 . This value is in slight contradiction 
with the recent observations of extragalactic Cepheids, which can be summa- 
rized as Ho > 60 kms~ 1 Mpc~ 1 . However, it is too early to speak about a 
serious conflict taking into account the many uncertainties and the physical 
possibilities (e.g. the Universe can be locally overdense - see the discussion in 
ref@). 

The age problem, if to take it seriously, can be solved with help of another 
successful cosmological model ACDM. This model assumes that SIq = 1 is 
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provided by the vacuum energy described by cosmological constant A and 
CDM. Using the limit on cosmological constant Qa < 0.7 and the age of 
Universe one obtains FIcdm > 0.3 and h < 0.7. Thus this model also predicts 
£lcDMh 2 ~ 0.15 with uncertainties 0.1. Finally, we shall mention that the 
CDM with f^o = ficDM = 0.3 and h — 0.8, which fits the observational data, 
also gives fth 2 « 0.2. Therefore Slti 1 ss 0.2 ±0.1 can be considered as the value 
common for most modelsd. 

We shall analyze here the candidates for CDM which naturally arise from 
elementary particle physics. The best known solution for strong CP violation 
implies axion, which can serve as CDM particle. The supersymmetrization of 
the theory, which includes axion, results in generalized neutralino - a linear 
combination of five neutral spin 1/2 particles (wino, bino, two higgsinos and 
axino, the fermionic partner of axion). This generalized neutralino is most 
natural candidate for CDM particle. And finally some attention will be given 
to the baryonic DM in connection with observations of machos. 

2 Machos and Baryonic Dark Matter 

The total number of microlcnsing events observed in the halo during last two 
years reached 10. Eight of them are observed by MACHO collaboration U 
and two - by EROS. The duration of lensing effect is determined by the lens 
mass. The distribution of observed durations yields the macho mass asB M = 
0.46^017^0- However, this value is model dependent. The most likely halo 
fraction of machos isQ / — 0.50io2o- The important result is observation El of 
45 microlensing events in Galactic bulge. For a given rotation curve the heavy 
bulge implies the lighter halo and thus the fraction of machos increases for a 
given number of observed events. 

The machos with these masses should be interpreted as white dwarfs. How- 
ever, to escape from the Hubble Deep Field Search these objects must be very 
faint, two magnitudes fainter than the disc white dwarf sequences. 

The DM in the halo of our Galaxy is found and most probably it is bary- 
onic. Could be all DM in the Universe only baryonic? This question is often 
asked nowdays (e.g. seeta). 

Let us discuss shortly the problems arising in the baryonic-dominated 
Universe. 

Nucleosynthesis requires fl™ cl < 0.02/i -2 . On the other hand the clusters 
provide > 0.2. Therefore, the baryonic density is small (fi^^ < S!^ M ), 

unless h < 0.3, which contradicts recent observations. If one arbitrary neglects 
this contradiction, the baryonic dominated Universe with f2s = ~ 0-2 

can be considered. Apart from IRAS data and POTENT analysis which give 
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S]w 1, such model faces serious cosmological problems, including the horizon 
and flatness problem and observed spectrum of fluctuations, which is impos- 
sible to explain without CDM and HDM. Probably, these problems could be 
solved in some artificial models with A term and vacuum defects (e.g. strings), 
but at present the corresponding calculations do not exist. 

As was mentioned above, the baryonic nature of machos have (or can have) 
the problems. .An interesting idea about the nature of machos was recently put 
forward in refLLj. These objects could be the neutralino stars, the formations 
produced by neutralinos and baryons around singularities EH in the distribu- 
tion of neutralino gas. The neutralino stars are produced naturally and they 
do not meet any problems connected with the Hubble telescope observations. 
Unfortunatelly, as was demonstrated in rcfO these objects produce too high 
gamma-ray flux due to annihilation of neutralinos. 

3 Axion 

The axion is generically a light pseudoscalar particle which gi¥££ natural and 
beautiful solutio n tp the CP violation in the strong interactionEil (for a review 
and references sec.il). Spontaneous breaking of the PQ-symmetry due to VEV 
of the scalar field <</>>= fpq results in the production of massless Goldstonc 
boson. Though fpQ is a free parameter, in practical applications it is assumed 
to be large, fpq ~ 10 10 — 10 12 GeV and therefore the PQ-phase transition 
occurs in very early Universe. At low temperature T ~ Aqcd ~ 0.1 GeV 
the chiral anomaly of QCD induces the mass of the Goldstone boson m a ~ 
Aqcd/ fpQ ■ This massive Goldstone particle is the axion. The interaction 
of axion is basically determined by the Yukawa interactions of field(s) <j> with 
fcrmions. Triangular anomaly, which provides the axion mass, results in the 
coupling of the axion with two photons. Thus, the basic for cosmology and 
astrophysics axion interactions are those with nuclcons, electrons and photons. 
Numerically, axion mass is given by 

m a = 1.9 • 10- 3 (iV/3)(10 10 GeV/fp Q ) eV, (2) 

where N is a color anomaly (number of quark doublets). 

All coupling constants of the axion are inversely proportional to fpQ and 
thus are determined by the axion mass. Therefore, the upper limits on emission 
of axions by stars result in upper limitsJpr the axion mass. In Table 1 we cite 
the upper limits on-axion mass from refEa, compared with revised limits, given 
recently by Raffeltlia. 

As one can see from the Table the strong upper limit, given in 1990 from 
red giants, is replaced by the weaker limit due to the horizontal-branch stars. 
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Tabic 1: Astrophysical upper limits on axion mass 



1990 



1996 



1 eV 



1 eV 



red giants 



1 • 10" 2 eV 



very 

uncertain 



hor. -branch 
stars 

SN 1987A 



not 

considered 
1 • 1CT 3 eV 



0.4 eV 
1 • 10~ 2 eV 



The upper limit from SN 1987A was reconsidered taking into account the 
nucleon spin fluctuation mN + N—*N + N + a axion emission. 

There are three known mechanisms of cosmological production of axions. 
They are (i)thermal production, (i) misalignment production and (iii) radiation 
from axionic strings. 

The relic density of thermally produced axions is about the same as for 
light neutrinos and thus for the mass of axion m a ~ 10~ 2 eV this component 
is not important as DM. 

The misalignment production is clearly explained in re fS 

At very low temperature T <C Aqqd the massive axion provides the mini- 
mum of the potential at value 9 = 0, which corresponds to conservation of CP. 
At very high temperatures T ^> Aqcd the axion is massless and the potential 
does not depend on 9. At these temperatures there is no reason for 9 to be 
zero: its values are different in various casually disconnected regions of the 
Universe. When T — > Aqcd the system tends to go to potential minimum 
(at 9 = 0) and as a result oscillates around this position. The energy of these 
coherent oscillations is the axion energy density in the Universe. From cosmo- 
logical point of view axions in this regime are. .equivalent to CDM. The energy 
density of this component is approximately E.M 



n a h 2 «2-(m a /10~ 5 eV)- 1A8 . (3) 

Uncertainties of the calculations can be estimated as 10 ±0 5 . 

Axions can be also produced by radiation of axionic strings^^. Axionic 
string is a one-dimension vacuum defect < 4>pq >= 0, i.e. a line of old 
vacuum embedded into the new one. The string network includes the long 
strings and closed loops which radiate axions due to oscillation. These were 
many uncertainties in the axion radiation by axionic strings (see refta for a 
review). Recently more detailed and accurate calculations were performed by 
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Battye and SchellardlLZI. They obtained for the density of axions 

n a h 2 « A{m a /10- 5 eV)- 1 - 18 (4) 

with A limited between 2.7 and 15.2 and with uncertainties of the order 10 . 
The overproduction condition fl a h 2 > 1 imposes lower limit on axion mass 
m a > 2.3 • 10~ 5 eV. Fig.l shows the density of axions f2 a /i 2 as a function 



Figure 1: Axion window 1996. The curves "therm." and "misalign." describe the thermal 
and misalignement productiwt* of axions, respectively. The dash-dotted curve corresponds 
to the calculations by Davis E3 for string production. The recent refined calculations EZI arc 
shown by two dashed lines for two extreme cases, respectively. The other explanations are 
given in the text. 

of the axion mass m a . The upper limits on axion mass from Table I are 
shown above the upper absciss (limits of 1990) and below lower absciss (limits 
of 1996). The overproduction region Vl a h 2 > 1 and the regions excluded by 
astrophysical observationsEj are shown as the dotted areas. 

The axion window of 1996 (shown as undotted region) became wider and 
moved to the right as compared with window 1990. The horizontal strip shows 
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Qcdm = 0.2 ± 0.1 as it was discussed in Introduction. One can see from Fig.l 
that string and misaligment mechanisms provide the axion density as required 
by cosmological CDM model, if axion mass is limited between 7 ■ 10~ 5 eV and 
7 • 10~ 4 eV. However, in the light of uncertainties, mostly in the calculations 
of axion production, one can expect that this "best calculated" window is 
between 3 ■ 10 -5 and 10 -3 eV. This region is partly overlapped with a possible 
direct search for the axion in nearest-future experiments (see Fig.l and refs£3). 

4 Generalized neutralino 

We shall consider now the supersymmetric theories where strong CP-violation 
is suppressed by PQ-symmetry (note, however, that in supersymmetric theories 
the new mechanisms for suppression of CP violation might appear, see e.g. E3). 
In supersymmetric theory the PQ symmetry breaking results in the production 
of the Goldstone chiral supermultiplet which contains two scalar fields and 
their fermionic partner - axino (a). The scalar fields enter the supermultiplet 
in the combination (fpQ + s) exp(a//pg), where s is a scalar field, saxino, 
which describes the oscillations of the initial field </> around its VEV value 
< <j) >= fpQ, and a is the axion field. This phase transition in the Universe 
occurs at temperature T ~ fpQ. As we saw in the previous section the axion 
is massless at this temperature and since supersymmetry is not broken yet, the 
axino and saxino are massless, too. The axion acquires the mass in the usual 
way due to chiral anomaly at T ~ Kqcd, while saxino and axino obtain the 
masses due to global supersymmetry breaking. 

The saxino is not of great interest for cosmology: it is heavy and it decays 
fast (mostly into two gluons). 

In the Minimal Supersymmetric Standard Model (MSSM) with broken 
PQ symmetry we have five spin 1/2 neutral particles: wino W3, bino B, two 
Higgsinos (Hi and H 2 ) and axino a. Generically the Lightest Supersymmetric 
Particle (LSP) is a linear superposition of these 5 fields: 

X = dW 3 + C 2 B + + C 4 H 2 + C 5 & (5) 

Further on we shall consider two extreme cases: generalized neutralino \ is 
strongly dominated by the axino state and it is dominated by the first four 
terms in eq. d). In the for mer case LSP is almost pure axino, in the latter- 
the usual neutralino. 

5 Axino 

How heavy the axino can be? The mass of axino is very model dependent. 
In the phenomenological approach, using the global supersymmetry breaking 
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parameter Msusy one typically obtains (e.g. 011) 

m~ a ~ M 2 SUSY /f PQ (6) 

For example, if global SUSY breaking occurs due to VEV of auxiliary field of 
the goldstino supermultiplet < F >= F g , then the axino mass appears due 
to interaction term (g j JpQ)aaF (F has a dimension M 2 ), and using < F >— 
F g = Mg USY one arrives at the value (||). ,— . 

The situation is different in supergravity. In refc2l the general analysis of 
the axino mass is given in the framework of local supersymmetry. It was found 
that generically the mass of axino in these theories is ma ~ ~ 100 GeV . 
Even in case when axino mass is small at tree level, the radiative corrections 
raise this mass to the value ~ m, 3 / 2 . This result holds for the most general form 
of superpotential.The global SUSY result, mj ~ 'm\/ 2 l fpQi can be reproduced 
in the local SUSY only if one of the superpotential coupling constants is very 
small, A < 10~ 4 , which implies fine-tuning. Thus, the axino is too heavy to be 
a CDM particle. _ 

The only exceptional case was found by Goto and Yamaguchi c3. They 
demonstrated that in case of no-scale superpotential the axino mass vanishes 
and the radiative corrections in some specific models can result in the ax- 
ino mass 10 — 100 keV, cosmologically interesting. This beautiful case gives 
essentially the main foundation for axino as CDM particle. 

The_cosmological production of axinos cap .occur through thermal pro- 
ductions or due to decays of the neutralinos c3e3. The axion chiral super- 
multiplet contains two particles which can be CDM particles, namely axion 
and axino. In this section we are interested in the case when axino gives 
the dominant contribution. In particular this can take place in the range 
2 • 10 9 GeV < fpQ < 2.7 • 10 10 GeV where axions are cosmologically unimpor- 
tant. 

Since axino interacts with matter vexy weakly, the decoupling temperature 
for the thermal production is very highcj: 

T d » 10 9 GeV(f PQ /10 u GeV). (7) 

Therefore, axinos are produced thermally at the reheating phase after inflation. 
The relic concentration of axinos can be easily evaluated for the reheating 
temperature Tr as 

"°- 6 ioofc^ ( f PQ } lo^u (8) 

Reheating temperature Tp < 10 9 GeV gives no problem with the gravitino 
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production. The relic density (||) provides SlcDAih 2 ~ 0.2 for a reasonable set 
of parameters ma, /pq and Tr. 

If the axino is LSP and the neutralino is the second lightest suneija^m- 
metric particle, the axinos can^ilso be produced by neutralino decays Ol3c3. 
According to estimates of reft3 the axinos are produced due to \ 5 + 7 
decays at the epoch with red-shift Zdec ~ 10 8 . Axinos are produced in these 
decays as ultrarelativistic particles and the free-streeming prevents the growth 
of fluctuations on the horizon scale and less. At red-shift z nr ~ 10 4 axinos 
become non-relativistic due to adiabatic expansion (red shift). From this mo- 
ment on the axinos behave as the usual CDM and the fluctuations on the scales 
A > (1 + z nr )ct nr (which correspond to a mass larger than 10 15 M Q ) grow as 
in the case of standard CDM. For smaller scales the fluctuations, as was ex- 
plained above, grow less than in CDM model. Therefore, as was observed in 
refE3, the axinos produced by neutralino decay behave like HDM. It means 
that axinos can provide generically both components, CDM and HDM, needed 
for description of observed spectrum of fluctuations. 

Unfortunatclly stable axino is unobservable. In case of very weak R-parity 
violation, decay of axinos can produce a diffuse X-ray radiation, with practi- 
cally no signature of the axino. 

6 Neutralino 

The generalized neutralino can be dominated by the first four terms in eq.(^|) 

X = C x W 3 + C 2 B + Czili + CiH 2 (9) 

i.e. by usual neutralino. 

The neutralino is a Majorana particle. With a unitary relation between 
the coefficients d the parameter space of neutralino states is described by 
three independent parameters, e.g. mass of wino M%, mixing parameter of two 
Higgsinos fj,, and the ratio of two vacuum expectation values tan/3 = f2 /vi- 
lli literature one can find two extreme approaches describing the neutralino 
as a DM particle. 

(i) Phenomenological approach. The allowed neutralino parameter space is 
restricted by the LEP and CDF data. In particular these data put a lower 
limit to the neutralino mass, m x > 20 GeV. In this approach only the usual 
GUT relation between gaugino masses, Mi : M 2 : M3 = a% : a 2 ■ 013, is used 
as an additional assumption, where cti are the gauge coupling constants. All 
other SUSY masses which are needed for the calculations are treated as free 
parameters, limited from below by accelerator data. 



9 



One can find Ihe relevant calculations within this approach in refs@E3 
and in the revie'wEj (see also the references therein). There are large ar- 
eas in neutralino parameter space where the neutralino relic density satis- 
fies flc'DAih 2 ~ 0.2 ± 0.1. This is especially true for heavy neutralinos with 
m x > 100 — 1000 GeV, refO. In these areas there are good prospects for 
indirect detection of neiitralinos, due to high energy neutrino radiation from 
Earth and Sun (seeE3S and references therein) as well as due to production 
of antiprotons and positrons in our Galaxy. The direct detection of neutrali- 
nos is possible too, though in more restricted parameter space areas of light 
neutralinos (see review E]). 

(ii) Strongly constrained models. This approach is based on the remark- 
able observation that in the minimal SUSY SU(5) model with fixed particle 
content, the three running coupling constants meet at one point corresponding 
to the GUT mass Mqut- Because of the fixed particle content of the model, 
its predictions are rigid and they strongly restrict the neutralino parameter 
space. This is especially true for the limits due to proton decay p — * K + v. 
As a result very little space is left for neutralino as DM particle. Normally 
neutralinos overdose the Universe (fl x > 1). The relic density decreases to the 
allowed values in very restricted areas where ^.-annihilation is accidentally 
large (e.g. due to the Z° exchange term - see refa. Thus, this approach looks 
rather pessimistic for neutraliniL as DM particle. 

In several recent worksE3~E3 less restricted SUSY models were considered 
with more optimistic conclusions about detection prospects. 

(in) Relaxed restrictions. In some recent works the restrictions described 
in (ii) are relaxed. In particular, ino the large number of models with relaxed 
conditions were analysed. It was found that for many models neutralian. can 
be discovered in the direct and indirect detection experiments. In refs. oEl the 
SUSY models only with basic restrictions were considered. 



7 SUSY models with basic restrictions 

Following refs.E3 ! i we shall consider here the restrictions to neutralino as DM 
particle, imposed by basic properties of SUSY theory. These restrictions are 
as follows: 

(i)Radiative Electroweak Symmetry Breaking (EWSB), which is considered 
as fundamental element of the analysis, (ii) No fine-tuning stronger than 1%, 
which is natural but very powerful requirement; it results in the upper limit 
to neutralino mass m x < 200 GeV, (iii) Restrictions from Renormalization 
Group Equations (RGE) and from particle phenomenology (accelerator limits 
on the calculated masses and the condition that neutralino is LSP), (iv) Limits 
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from b — y sj decay taken with the uncertainties in the calculations of the decay 
rate and (v) 0.01 < fi x /i 2 < 1 as the allowed relic density for neutralinos. 

At the same time some restrictions are lifted as being too model-dependent: 
(i) No restrictions are imposed due to p —y Kv decay, (ii) Unification of cou- 
pling constants at the GUT point is allowed to be not exact (it is assumed 
that new very heavy particles can restore the unification), (iii) unification in 
the soft breaking terms is relaxed. Following refill it is assumed that masses 
of Higgses at the GUT scale can deviate from the universal value mo as 

m 2 Hi {GUT) = m 2 (l + Si) (i = l,2). (10) 

This non- universality affects rather strongly the properties of neutralino as 
DM particle: the allowed parameter space regions become larger and neutralino 
is allowed to be Higgsino-dominated, which is favorable for detection. 

Some results obtained in refe^ta are illustrated by Figs. 2-3. 

In Fig. 2 the regions excluded by the LEP and CDF data are shown by dots 
and labelled as LEP. The regions labelled " fine tuning" have an accidental com- 
pensation stronger than 1% and thus are excluded. No-fine-tuning region in- 
side the broken-line box corresponds to a neutralino mass m x < 200 GeV. The 
region "EWSB+particle phenom." is excluded by the EWSB condition com- 
bined with particle phenomenology (neutralino as LSP, limits on the masses of 
SUSY particles etc). In the region marked by rarefied dotted lines neutralinos 
overdose the Universe (il x h 2 > 1). The solid line corresponds to Too = 0. The 
regions allowed for neutralino as CDM particle (0.01 < fl x h 2 < 1) are shown 
by small boxes. As one can see in most regions the neutralinos are overpro- 
duced. The allowed regions correspond to large \X annihilation cross-section 
(e.g. due to Z°-pole). 

Fig. 2a and Fig. 2b differ only by universality: in Fig. 2a 8\ = 82 = (mass- 
unification), while in Fig. 2b Si — —0.2 and 82 = 0.4. The allowed region in 
Fig. 2b becomes much larger and is shifted into the Higgsino dominated region. 
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Figure 2: The neutralino parameter space for (a) mass-unification case Si = 62 = and 
(b) for non— universal case <5i = —0.2, &2 = 0.4. Both cases are given for tan (3 = 8. 

Let us discuss now the predictions of this model for direct and indirect 
detection of neutralinos. Direct detection is based on observarions of recoil 
nuclei from neutralino-nucleus scattering. As indirect detection we shall con- 
sider here the registration of high energy neutrinos from neutralino- neutralino 
annihilation in the center of Earth and Sun. 

In Fig. 3a the scatter plot for the rate of direct detection with the Ge 
detectorcj is given for the non-universal case (<5i = 0, 62 = —0.2) and tanJi = 
53. We notice that, for some configurations, the experimental sensitivity l3 is 
already at the level of the predicted rate. 

In Fig. 3b we show predictions for the updoing flux of muons produced by 
neutrinos from neutralino-neutralino annihilation in the core of Earth. The 
muon flux from the direction of the Sun is shown in Fig. 3c. The. horizontal 
solid curves in both cases present the observational upper limits E-l One can 
see that in examples given above the muon fluxes can be reliably detected by 
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future gigantic neutrino telescopes. 



Figure 3: Direct (a) and indirect (b and c) detection of neutralinos for tan/3 = 53 and 
different ncutralino masses. Fig. 3a presents counting rate in units events / (kg ■ day). Fig. 3b 
gives the underground muon flux from the Earth-core direction for <$i = 82 = (solid line), 
81 =0,82 = —0.3 (dashed line), and <5i = 0.7, 82 = 0.4 (dotted line). Fig. 3c gives the 
underground muon flux from the direction of the Sun for the case 81 = and £2 = —0.3. ; 
the fluxes are given in units cm~ 2 s~ 1 . 



8 Conclusions 

The baryonic DM is discovered in the halo of our Galaxy. Machos have a mass 
between 0.1 — 1.2M© and comprise between 10 — 100% of the total mass of 
galactic halo. However, the bulk of DM observed in the Universe can hardly 
be dominated by baryons. 

The most successful cosmological models require CDM with density Vth 2 « 
0.2 ±0.1. The minimal supersymmetric extension of SM, with strong CP vi- 
olation suppressed by PQ-symmetry, predicts the generalized ncutralino as 
a superposition of usual neutralino and axino. The Goldstone chiral super- 
multiplet contains axion and axino. Therefore, in this model there are three 
natural candidates for CDM particle: axion, axino and neutralino (or a linear 
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superposition of the latter two). 

The new axion window corresponds to axion masses between 3 • 10~ 5 and 
3 • 10~ 3 eV, i.e. it only partly overlaps with the range of search in microwave 
cavity experiments. 

Axino can provide both CDM and HDM. The direct observation of this 
particle seems to be impossible. 

Neutralino remains most attractive CDM candidate. In the models with 
radiative EW symmetry breaking the properties of neutralino are restricted, 
but there are many configurations where the neutralino can provide the re- 
quired flcDAih 2 and can be found by direct and indirect methods. 
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